Dynamics of quantum Hall stripes in double-quantum-well systems 



R. Cote 1 , H.A. Fertig 2 
1 Departement de Physique and CERPEMA, Universite de Sherbrooke, 
Sherbrooke, Quebec, Canada J1K 2R1 
2 Department of Physics and Astronomy, University of Kentucky, 
Lexington KY 40506-0055 
(February 6, 2008) 

The collective modes of stripes in double layer quantum Hall systems are computed using the 
time-dependent Hartree-Fock approximation. It is found that, when the system possesses spon- 
taneous interlayer coherence, there are two gapless modes, one a phonon associated with broken 
translational invariance, the other a pseudospin-wave associated with a broken (7(1) symmetry. For 
large layer separations the modes disperse weakly for wavevectors perpendicular to the stripe ori- 
entation, indicating the system becomes akin to an array of weakly coupled one-dimensional XY 
systems. At higher wavevectors the collective modes develop a roton minimum associated with a 
transition out of the coherent state with further increasing layer separation. A spin wave model of 
the system is developed, and it is shown that the collective modes may be described as those of a 
system with helimagnetic ordering. 
PACS: 73.20.Mf,73.21.Fg,71.10.Hf 



I. INTRODUCTION AND SUMMARY OF RESULTS 

Quantum Hall systems are by now well-known to support a broad range of condensed matter phenomena jjj . These 
include localization physics, electron solidification (Wigner crystals), chiral Luttinger liquids at quantum Hall edges, 
and Fermi liquid behavior. Recently added to this list is highly anisotropic transport in moderate magnetic fields, in 
which electrons populate several Landau levels Such behavior was anticipated by mean-field studies predicting 
unidirectional charge density wave (CDW) ground states for electrons in high Landau levels While more 

careful studies || demonstrate that within Hartree-Fock theory such ground states are unstable to the formation of 
"modulated" stripe states which are essentially highly anisotropic two-dimensional Wigner crystals ( "stripe crystals" 
0), quantum fluctuations may restore the translational symmetry along the stripes ]7|-|l0|]. It has been noted that 
the CDW state has the symmetry of a two-dimensional smectic M, and this analogy with liquid crystals has been 
exploited to yield a number of useful results |ll[ . 

A seemingly unrelated set of physical phenomena occurs in double quantum well systems (DQWS) in the quantum 
Hall regime. For large magnetic fields, these systems support spontaneous interlayer coherence [fl2||, and exhibit an 
associated Goldstone mode [[l3 
these systems, contributing to recent excitement about them. 

When immersed in more moderate magnetic fields, such that several Landau levels are occupied, it becomes evident 
that the physics of stripes and interlayer coherence are simultaneously relevant in the bilayer system. This situation 
was recently explored |Iq ]. For a system in which each layer is half-filled in the Nth Landau level, with some 
simplifying assumptions (lower Landau levels filled and inert, Landau level mixing negligible, and Zeeman coupling 
large enough that spins are polarized) it was found that: (1) For small enough layer separation, and any value of N, 
an interlayer coherent state of uniform density is formed. (2) Above a critical layer separation d c (N), a unidirectional 
density modulation sets in such that the electron density oscillates between the layers as a function of position Jl9| . 
This state is simultaneously interlayer coherent and smectic in spatial symmetry. With increasing d, the regions of 
coherence become relatively narrow in one direction, forming "linear coherent regions" (LCR's). (3) For negligibly 
weak interlayer tunneling, the coherence is spontaneous and a corresponding Goldstone mode should be present in the 
excitation spectrum. When the LCR's are sufficiently narrow compared to their separation (which occurs at large d 
and/or large N), the system may be thought of as an array of coupled one-dimensional XY models. Different possible 
quantum-disordered states may exist [ [18[|2C(| , which in principle can be distinguished by tunneling experiments. 

In order to predict which, if any, of these quantum disorderings might occur, it is necessary to understand the 
low-energy collective modes in some detail, and to create a model that captures them which is sufficiently simple to be 
susceptible to further analysis Q . This motivates our present study. Beyond its use in specifying a low-energy theory, 
the collective modes are interesting on their own, as they are sometimes detectible in electromagnetic absorption |2l| l 
or Raman scattering E^] experiments, and allow one to learn about which states are realized in real samples. 

In what follows, we compute the collective mode spectrum of the electrons in a DQWS, in a perpendicular magnetic 
field, in which several Landau levels are occupied, using the time-dependent Hartree-Fock approximation (TDHFA) 



14 . A precursor of the Josephson tunneling |Tq,M may have been observed JT71] 
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p3| . Counting both the spin and layer index degrees of freedom, the situations we consider involve AN filled Landau 
levels, which are taken to be inert (as appropriate for large cyclotron energy huj c , where lu c = eB/m*c, B is the 
magnetic field, and m* the effective mass of electrons in the quantum wells, and large Zeeman coupling.) For 
simplicity we ignore finite thickness of the wells. In addition to the AN filled levels, there is a level, with Landau index 
N, which contains precisely enough electrons so that the total filling factor v, defined as the number of electrons per 
magnetic flux quantum passing through the plane, is AN + 1. The system considered has no interlayer bias, so that 
the highest level in each layer will be half-filled on average. 

As discussed above, if the individual layers are isolated, this is precisely the situation in which one expects stripes 
to form. For finite layer separations, it is clear that the ground state will have the electron-rich regions in one layer 
align directly above the hole-rich regions in the other. Furthermore, for a range of layer separations, the stripe edges 
hybridize between the layers so that the electron occupation continuously shifts between the layers as a function of 
position Jig . An example of such a ground state is illustrated in Fig. ||. To understand this state and its collective 
modes it is convenient to adopt a pseudospin representation, in which an electron pseudospin is "up" when occupying 
one layer, and "down" when occupying the other. In this language, the stripe ground state involves a unidirectional 
spatial tumbling of the spin, and may be thought of as having helimagnetic ordering (2J]. As the separation increases, 
the electron density tends to reside more sharply in one well or the other, and the transition regions between layers 
become narrower |18| . 

In the absence of interlayer tunneling, the pseudospin density may tumble from +z to —z by passing through the 
x direction, the y direction, or any direction in between. Thus, each of the transition regions - the LCR's mentioned 
above - supports an XY degree of freedom, which is spontaneously broken in the ground state. In addition, this state 
obviously has broken translational symmetry. 

Because of the two broken symmetries, there are two gapless excitations above the stripe ground state. One 
is associated with the broken XY symmetry - alternatively thought of as a phase coherence between the wells 
[ p2| , ^3| - which disperses linearly with wavevector. The other is a phonon-like mode which disperses quadratically 
for wavevectors parallel to the stripes, and linearly perpendicular to them. This is analogous to what is seen in 
incommensurate helimagnets p4j . 

A typical dispersion relation for the XY mode is displayed in Fig. || for several values of d/£, where £ = y/ hc/eB 
is the magnetic length. As may be seen, the linear mode tends to disperse more strongly along the stripes than 
perpendicular to them, and in the limit of large d, the latter dispersion may become quite weak. This may be 
understood in terms of the exchange coupling among LCR's, which vanishes as they become arbitrarily narrow |L8| . 
It is possible that quantum fluctuations can effectively wipe out the inter-LCR exchange coupling, leading to a state 
analogous to a "sliding XY model" Jl8| , p5[ . 

Several other features appear in the collective mode spectrum. As is typical of such calculations, a number of higher 
energy modes are present, one of which may be interpreted as a gapped "out-of-phase" phonon mode, in which the 
stripes in different layers oscillate against one another. In Fig. we see that this mode becomes degenerate with 
the XY mode in its dispersion relation along the stripe direction and develops a roton minimum which touches zero 
at a critical separation dM, signaling an instability in which modulations form along the stripes. The instability, 
however, is first-order in nature since the modulation amplitude changes discontinuously at the transition |l8| . It is a 
surprising feature of our results that a collective mode appears to go soft so close to this transition, and this behavior 
in principle allows a detection of the transition by inelastic light scattering. Whether such a bilayer stripe crystal is 
stable with respect to quantum fluctuations analogous to those being discussed in the single layer case ]7|-^0|] is at 
present unclear. 

This article is organized as follows. In Section II we sketch the Hartree-Fock and TDHFA formalism, providing 
some details about its application to the two-layer stripe system. Section III contains a more detailed description of 
our results, and Section IV shows how the low-energy spectrum may be understood in the pseudospin language. We 
conclude with a summary and discussion in Section V. 



II. HARTREE-FOCK DESCRIPTION OF THE CDW GROUND STATE 



We consider an unbiaised symmetric DQWS in a perpendicular magnetic field B = Bz at total filling factor 
v = AN + vq where N = 0, 1, 2, ... is the Landau level index and vq is the filling factor of the partially filled level. 
Each Landau level has two spin states and two layer states specified by the index j = R,L. The layer states hybridize 
into symmetric (S) and antisymmetric (AS) states in the presence of tunneling. We assume the magnetic field to be 
strong enough so that the lower AN levels are completely filled with electrons and can be considered as inert. The 
Zeeman energy is assumed to be much larger than the S-AS gap and so there is no spin texture in the ground state. 
The electron gas is completely spin polarised and only the layer degree of freedom need be considered. In the ground 
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state, the electronic charge is equally distributed between the two wells. 

In the Landau gauge where the vector potential A = (0, Bx, 0), the electron wavefunctions are given by 

Viv.x^r) = JL e -«»/<V JV (x - X) Xj (z), (1) 



V L v 

where N and X are the Landau level and guiding center indices and Xj{ z ) lSi the envelope wave function of the lowest- 
energy electric subband centered on the right or left well. ipN {x) is an eigenfunction of the one-dimensional harmonic 
oscillator. Because the fully occupied Landau levels are considered as inert, we need only consider the partially filled 
Landau level. We can then drop the Landau level index from now on. 

To describe the various charge density wave (CDW) ground states, we define the operator 

= ±J2 e - iqMqye2/2 ix^x-^, (2) 
"<p x 

where i, j = R,L and N v is the Landau level degeneracy. In the ground state, (pj,j(q)) ^ where q is a reciprocal 
lattice vector of the CDW. By definition {pr,r(<\ = 0)) = (pl,l{<1 = 0)) = vo/2. In the Hartree-Fock approximation 
(HFA), the ground-state energy per electron in units of e 2 /n£ in the partially filled level can be written in term of 
these operators as 

Ehf = ~~f Re [{pr,l (q = 0))] (3) 



+^ E \ H (q) - x \\{pr,r (q))l 2 + \{PL,L (q))i 2 



q 



+ 77J2 H ^ (Pn,n(ci)) ( P l,l (-ci)) 



-Yx(q) (p fl , L (-q)><p L ,fl(q)>. 



^ q 

In this equation, f is the tunneling energy in units of e 2 /k£, and i?(q) ,X (q) are the Hartree and Fock infra-well 
interactions while H (q) and X (q) are the Hartree and Fock infer-well interactions. For very narrow wells where 
Xj(z) are highly localized, these interactions are given by 

^(q) = (ji) [L° N (?e/2)] 2 e-M*, 

H(<l) = (£j [L%( q H>/2)] 2 e-^/>e-« d , (4) 

X(q) = S™dyJ {yqt)e-y 2 /->[L%(y*/2)] 2 , 

X(q) = j °° dyJ (y q £)e-v 2 /\-y^ [L% ( y y2)] 2 , 

where d is the center-to-center separation between the wells and LP N (x) are generalized Laguerre polynomials. Because 
of the neutralizing positive backgrounds of ionized donors on both sides of the DQWS, we have H (0) = H (0) = in 
Ehf- 

The order parameters {(pij (q))} are computed by solving the HF equations of motion for the single-particle Green's 
function 



dj (X, X 1 , t) = - (T Ci , x (r) c\ x , (0)) (5) 

whose Fourier transform we define as 

G id (q,r) = -L]T e- L M^')5 x . x ,_ qyl ,G t , 1 (X,X',r) , (6) 

4 X,X' 

so that dj (q,r = 0~) = {p jti (q)) . 

In the HFA, these equations of motion are given by 

( ( ifilJ n + A 4 ) Sq.q" ~ U R,R (q, q") ^q.q" ~ U R.L (q, q") \ m 

^ \ tfq.q" - U l ,r (q,q") (itkJ n + p) 5 q , q „ - C/ L , L (q,q") ) ( '> 



q" 



/ Gr.r (q", w n ) G^.l (q", u> n )\_{ M q . 
^Gi, fl (q",a; n ) G L , L (q'>„) J \ M q , 
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where p is the chemical potential. The effective potentials Uij (in units of e 2 / k£) are defined by 

Ur,r (q, q') = [H (q - q') - X (q - q')] ( P r,r (q - q')> 7q,q' + # (q - q') (Pl,l (q - q')) 7q,q' , (8) 

(q, q') = \H (q - q') - X (q - q')] (p L , L (q - q')> 7q,q' + H (q - q') ( P r,r (q - q')> 7q,q' - (9) 

(q, q') = -X (q - q') (p i>fl (q - q')) 7q,q' , (10) 

U L ,R (q, q') = -X (q - q') {pr,l (q - q')> 7q,q', (H) 

where 7 q , q ' = e~ lqxq '^/ 2 . The procedure to solve for the (pij (q))' s is described in detail in Ref. |23|j . 

It is instructive at this point to describe the electron state in a DQWS by using a pseudospin language where states 
R and L are mapped to up and down spin states. For this, we define the pseudospin operators 

(12) 



P(q) 


i 

~ 2 




(q) + PL,L (q)] , 


Px (q) 


1 

~ 2 


[pR,L 


(q) + PL,R (q)] , 


Py (q) 


1 

~ 2i 


[PR,L 


(q) - PL,fl (q)] 


Pz (q) 


1 

~ 2 


[pR,R 


(q) - PL,L (q)] , 


p± (q) 


= Px 


(q)x 


+ Py (q) y- 



These operators act in the restricted Hilbert space of the partially filled Landau level. To get a real space representation 
of the ordered states, we take the Fourier transform 

(Pij (r)) = E 0*J (<!)) elq ' r - ( 13 ) 
q 

A local filling factor in each well can then be defined as 

Ui (r) = <p M (r)) . (14) 

For unidirectional modulations, the quantity 2/j (r) will take values between and 1. Similarly, a Fourier transform 
of the order parameters defined in Eq. ( [l^ ) will define the total density (p(r)) and spin density (p M (r)). Note that 
these densities are related to the density of orbit centers and not to the real density of electrons. For instance, when 
Landau mixing is neglected, the real "densities" n^j (q) are given by 

n id (q) = N^Fn (q) Pilj (q) (15) 

where (q) is a form factor appropriate to Landau level TV given by 

F N (q) = exp ( d?L.) L% ( f?) . (16) 



4 J " \ 2 

In the pseudospin language, the HF energy can be rewritten as 



Ehf = — — i Re [{ Px (0)}] (17) 



+-E T W (p(-q))(p(q)> 

+— yy.(q) <p.(-q))<p* (q)> 

q 

+-£^(q) W(-q))-(p± (q))- 



^0 

q 
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where the effective interactions are defined by 



T(q) = H (q) + H (q) - X (q) 
J z (q) =if^(q)- J ff(q)-X(q) 
Jj_ (q) = -X (q) . 



(18) 



We derive the dispersion relations of the collective excitations of the CDW states in the DQWS by tracking the 
poles of the retarded density and pseudospin response functions. These are obtained by analytical continuation of the 
two-particle Matsubara Green's functions 



Xi,j,k,t (q, q'; t) = -N+ (Tp id (q,r) p k j (-q', 0)} + {p id (q)> (p kJ (-q')) 



(19) 



which we compute in the generalized random phase approximation (GRPA). The procedure is explained in details in 
Refs. |fii|]23f|. It is convenient to work in the pseudospin language where we can define the matrix 



/ Xp,p Xp,x Xp.y Xp,z 

Xx.p Xx.x Xx,y Xx.z 

Xy-p Xy,x Xy,y Xy,z 

\ Xz,p Xz,x Xz,y Xz,z 



(20) 



These pseudospin Green's functions are related to the original Green's functions of eq. (|T^) by the transformation 

(21) 



i '. 1 1 I' i' 



where v, p 



p, x, y, z. ct,- • are Pauli matrices and a p — 



1 
1 



The summation of the bubbles and ladder diagrams of the GRPA can be expressed as an equation of motion for \ 
of the form 



(w + iS) * (q, q', a;) - £ F (q, q") X (q", q» = D (q, q') . 



(22) 



Here F and D may be written schematically as 



F 



-2i 



(p) sin a [T — T'] (p x ) sina [J± — J'j] (p y )sma[J± 

(p x ) sina [Jj_ — T'] (p) sina [T — J' ± ] (p z )cosa[J z 

(p y ) sina [Jj_ — T'] — (p z ) cosa [J z — J\] (p) sina [T - 

(p z ) sina [J z - T'] (p y ) cosa [Jj_ - J' ± ] t - (p x ) cos a [J± - J' ± ] 



-J'±. 
J',1 



(p z ) sina [J z - J' z ] 
— (p y ) cosa [Jj_ — J^] 

* + (Px) cos a [J± " J' z ) 
(p) sina[T - J' z ] 



and 



i(p)(sina) i{p x }(sina) 

i (Px) (sina) i(p)(sina) 

i (p v ) (sina) —i(p z )cosa 

i(p z )(sma) i (py) (cosa) 



i (p y ) (sina) 
i (p z } cos a 
« (p) (sin a) 

-i (p x ) (cosa) 



i (p z ) (sina) 
— i (p y ) (cosa) 
i (p x ) (cosa) 
i (p) (sin a) 



(23) 



(24) 



These matrices contain wavevector dependence that is not shown. For example, in the F matrix, the first entry is 
explicitly given by 



2i (p) sin a [T - T'] -> -2i (p (q - q")) sin 



q x q"Z 



inl n 



[T (q-q")-T(q") 



(25) 



The other terms in the matrix, as well as those in the matrix D, have analogous definitions. Finally, for the tunneling 
term in F 34 and F43 is diagonal in the wavevector index i.e. t — > tSq^". Eq. (|2^) may be solved by diagonalizing F 
p3| , and the collective mode frequencies found from its eigenvalues. From the eigenvectors of the F matrix, it is also 
possible to extract the motion of the guiding-center densities and of the pseudospin in a given mode. 

It is interesting to note that for (p x . y ) — and t — 0, the forms of F and D indicate that in the equations of motion 
Xp,p> Xp.z, Xz, P , Xz,z completely decouple from Xxx, Xxy, Xyx, Xyy This indicates that distortions of the stripes 
involving motion of charge cither within the layers or between them is completely decoupled from any "in-plane" XY 
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motion of the pseudospins; i.e., phonon modes and spin-wave modes will create poles in different, distinct response 
functions. The presence of coherence - a non-vanishing (p x ) or (p y ) - mixes these motions, so that poles from all the 
collective modes appear in all the response functions. This phenomenon is closely related to "spin-charge coupling" 
that is generically present in multicomponent quantum Hall systems Q, and has important consequences for the 
charged excitations in this system |l8|. For the collective modes, we will see below that the low-energy interlayer 
charge degrees of freedom (distortions that change p z ) are in a sense conjugate to the in-plane degrees of freedom 
{px,y), and distortions of both are involved in any given collective mode. 



III. NUMERICAL RESULTS 

In Ref. jl8| the energies of several ordered ground states at vq = 1 were computed in the HFA. The states considered 
there were a uniform coherent state (UCS), a unidirectional coherent charge density wave (UCCDW), a modulated 
unidirectional charge density wave (MUCDW) and a coherent Wigner crystal (CWC) with a square lattice. We refer 
the reader to Ref. for more detailed discussions of these states. For N > 0, it was found that the ground state of 
the electron gas evolves from the UCS for d < d c (N) to the UCCDW at larger values of d and finally to the MUCDW 
as the separation between the wells increases. At large d, the Wigner crystal state is only lowest in energy for Landau 
level N = 0. (There is however a small region of d for d > d c in N = for which the UCCDW is lowest in energy (see 
Ref. |l8|])). For all ordered states, the lowest energy is obtained when the density pattern in both wells are shifted 
with respect to one another. Moreover, coherent states (states with non-zero value of (pr^l (q)))j when they exist, 
have lower energy than their incoherent counterparts. For the MUCDW it was impossible to find a coherent version 
in the region where it has lower energy than the UCCDW. An interesting result of the HF calculation is that the 
local coherence (pr,l (•*)) is maximal when the charge density is equally shared by both wells. For the UCCDW, this 
occurs along channels called linear coherent regions (LCR). As the separation between the wells increases, the width 
of the LCR's becomes very small. Figs. 1-3 summarizes the HF results for level N = 2. Fig. [l] shows the energy of 
the four states defined above as a function of d in Landau level N = 2 and in the absence of tunneling. (We remark 
that all results presented in the present paper are for vq = 1. Also, in the absence of tunneling, the phase of (pr,l (r)) 
is arbitrarily choosen so that all spins point in the x direction in the ground state). As a measure of the coherence 
of a given state, we use (pr,l (0)} . This quantity takes its maximal value 1/2 (at Vq = 1) in the UCS. In Fig. |l|, the 
coherence decreases slowly for the UCCDW but very rapidly for the CWC and is essentially zero in the MUCDW. 

In Landau level N = 2, the UCS is lower in energy for d/£ < 0.7 at which point it evolves continuously into the 
UCCDW. At djl ~ 1.6, there is a first order transition into the MUCDW. Fig. || shows the guiding- center density 
in the right and left wells as defined in Eq. ( |l4| ) and the pseudospin pattern (Eq. ([l2"|)) for the UCCDW at djl = 0.7 
and djl — 2.0. The formation of the LCR's is clearly visible in this figure. At large d, the coherence is very small 
(see Fig. and the guiding-center densities approach the stripe pattern appropriate to decoupled layers with filling 
factor 1/2. For completeness, we also show in Fig. || the evolution of band structure jl^] E (X) and of the density of 
states in the UCCDW at N = 2 and for t = 0. 

We now consider the collective excitations of the UCS and UCCDW. The dispersion is obtained by solving Eq. (|2S 
for the suceptibility \- I n the UCS, we can easily solve this equation to get the response functions 



/0 





\0 



-4aJ> 

— 1 

-4af> 



j*-4ab 
b 



-4ab 



(26) 



where 



(Px (0)) 



(27) 



and 



a (q) = t - ( Px (0)) [Jl (0) - J z (q)] , 
&(q) =*-<A e (0)>[Jl(0)-J±(q)]. 



(28) 
(29) 



The dispersion relation of the collective mode of the UCS is given by 



2 V a (q)b(q). 



(30) 



This mode is a Goldstone mode (at t = 0) associated with the broken XY symmetry of the UCS. For small wavevectors, 
it disperses linearly in q for t — and d ^ 0. This coherence mode represents an elliptical motion of the pseudospins 
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around the x axis. The pseudospin motion becomes more and more confined to the x — y plane as g decreases. For 
level N = 2 and in the absence of tunneling, the dispersion relation of this mode softens at interlayer separation 
d/l f» 0.64. This softening occurs at a finite value q£ = 0.92 of the wavevector signaling the onset of the formation of 
the UCCDW state with a wavelength (separation between stripes in a given layer) of approximately X/l — 27r/0.92. 

For our choice of phase, the pseudospins in the UCCDW rotate in the x — z plane. This implies that (p y (q)) = 0. 
Moreover, in this shifted state, there is no modulation of the total density, so that (p(q)) = as well. The F and D 
matrix introduced in Eqs. (E3L) and (El) then simplify to 



F = -2i 



(p x ) sina [J± - J' ± ] (p z ) sina [J z - J' z ] 

(p x ) sin a [J_l - T'] (p z ) cosa [J z - J' ± ] 

- (p g ) cos a [ J z - J' x ] t + (p x ) cosa [J± - J' z ] . 

(p z ) sina [J z - T'] t - (p x ) cos a [J± - J' ± ] / 1 ' 



and 



i (p x ) (sina) i (p z ) (sina) 

U I i{Px) (sina) i(p z )cosa , . 

' —i(p z )cosa i(p x )(cosa) ' ' "~ ' 

i(p 2 )(sina) —i(p x )(cosa) 

As is clear from Eq. ([H]), The interlayer coherence present in the UCCDW introduces a coupling between the lon- 
gitudinal and transverse response functions. From Eqs. ( |3l"| ) and (|32|), we see that x pp — > as fc| — » so that 
the coupling with the density response gets very small for small wavevector parallel to the stripes in which case the 
response is dominated by the pseudospin motion. Fig. [| shows the dispersion relations of the lowest four collective 
modes of the UCCDW at d/l — 1.0 with and without tunneling. The dispersion is given for wavevector along the 
direction of the stripes with k± = 0. These curves are obtained by tracking the poles in the four response functions 
Xp,pi Xx,xi Xy,yy Xz,z for wavevector k along the desired direction in the Brillouin zone. From the weight of given pole 
in each response function, we can infer the nature of the mode. The low-energy dispersion consists of an in-phase 
and out-of-phase phonon modes (empty squares in Fig. ^]) that both involve a coupling between the density p and 
pseudospin p x . The in-phase phonon mode is gapless while the out-of-phase phonon is gapped. Both phonons are 
gapped for fc|| = 0, \kj_\ > (see Fig. ^|(b)) in contrast with what happens for stripes in single quantum well sys- 
tems where the phonon frequency vanishes for all fcii = 0. These behaviors are distinct because the nature of the 
interstripe-coupling in the single layer and double layer systems is different in an important way. In the single layer 
system, there is very little exchange interaction between stripes. Any dispersion in the perpendicular direction comes 
from direct coupling, i.e., the Hartrcc interaction. In the single layer case, modulations along the stripes are present 
and in principle introduce a gap in this direction. In practice, the energy cost for "sliding" stripes with respect to 
one another is small because the modulations are weak, and is nearly averaged out due to the long-range nature of 
the Coulomb potential. Thus, in calculations such a gap is essentially immeasurable g. In the present case, the 
coupling between stripes is due to exchange, it is present even in the absence of modulations along the stripes, and is 
not averaged away due to the long-range nature of the interaction. The exchange coupling is set by matrix elements 
between single-particle states in different LCR's ; these become small in the limit of large d or N but in general 
are not negligible. 

In addition, one may also clearly see a linearly dispersing gapless XY mode in Fig. |^ which, as in the UCS, 
represents a motion of the spins in the y — z plane. In the UCCDW, the dispersion relation of this mode is folded 
into the first Brillouin zone (as is the case for the phonon modes as well). In Fig. [|, we show two branches of this 
XY mode represented by the filled diamonds. Fig. ^(b) shows how tunneling affects these dispersion relations. As 
expected, the phonon modes are not dramatically affected by switching on the tunneling while the phase (XY) mode 
becomes gapped. 

The dispersion relations of the phase and in-phase phonon modes are plotted in Figs. || and ^ for directions parallel 
and perpendicular to the stripes and for several values of d/l. For the phase mode, the dispersion is linear in both 
directions but weaker in the perpendicular direction. Comparing Figs. ^|(a) and (b), we see that the phonon dispersion 
is quadratic along the stripes and linear for direction perpendicular to the stripes. As d increases, the phonon and 
phase mode dispersions in the perpendicular direction become very weak and eventually their gaps vanish in the 
limit of very large d. The suppression of these gaps reflects the shrinking of the exchange coupling discussed above, 
and indicates that the system is essentially an array of weakly coupled one-dimensional systems in this limit. For 
the phonon mode, our results are consistent with the calculated dispersion for the phonon mode of stripes in single 
quantum well Q.The qualitative behavior of the gapless modes at low energies may be understood in terms of a 
spin-wave model which will be developed in the next section. 
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In Fig. ||(a), the out-of-phase phonon mode is seen to become degenerate with the XY mode at large values of k\\. 
Both modes soften at approximately q»/ f^ZLJ = 3.1 when d increases and at d/£ — 1.6 (a is the separation between the 
stripes in a given layer) they become unstable. The period of modulation along the stripes implied by this instability 
is consistent with the formation of a MUCDW or highly anisotropic Wigner crystal with one electron per unit cell in 
each well. The softening apparently accompanies a first order transition into the MUCDW, since both the UCCDW 
and MUCDW exist as solutions to the HFA both above and below the critical separation, and cross in energy very 
close to it. Note that for stripes in a single quantum well system, within HF theory the unidirectional CDW state is 
always unstable with respect to the formation of an anisotropic Wigner crystal. Here, the instability only occurs at 
large enough values of d. In the MUCDW, the coherence is quickly lost with increasing d as can be seen on Fig. [jj 



IV. SPIN WAVE ANALYSIS 



As mentioned previously, the symmetry of the ground state and the low-energy excitations are formally quite similar 
to those of a non-collinear ferromagnet, with helimagnet ordering ]24| ]. In this section, we demonstrate that such a 
model can be constructed that captures the low-energy behavior, and explicitly demonstrates the origin of the two 
low-energy modes. Our analysis uses the pseudospin analogy introduced above but now with the real electron density 
difference between wells which we denote by 



Sz(t) 



tpR ( r ) ( r ) - V4, (r) V>l (r) 



With the natural definitions for spin raising and lowering operators, 

S+(r) = S- (r) f =^ (r)V-L (r) 



(33) 



(34) 



we then have in-plane spin components S x = ^ [S+ + S-], S y = A [S+ — S-]. These spin operators obey the usual 
commutation relations 



[S i (r),S j (r')] = y2^kS k (r)S(r-r') 



k 



(35) 



where k = x,y, z and e%jk is the antisymmetric tensor. 

These spin operators are obviously related to the spin operators defined in Section II. The connection is most easily 
seen when the spin commutation relations are Fourier transformed, to give 

[s,(q),Mq')] = ^(q + q')- (36) 

This should be compared to the guiding center density and spin density operators [Eqs. (fL2|)] which obey the algebra 

N<f, [p (q) , p (q')] = ~i sin (q X q'l 2 /2) p(q + q') 
[p (q) , p^ (q')] = ~i sin (q X q'l 2 /2) Pi (q + q') , 



[Pi (q) . Pi (q')l = i£ijk cos (q x q'l 2 /2) Pk (q + q') ■ ^ 3 



(37) 



It is clear that in the limit of small q and q' the density and spin density operators decouple. Moreover, if we make the 
identification Si (q) ~ N^pi (q) one can see the direct connection between the microscopic operators and the effective 
ones used in this section in the long- wavelength limit [ p6[ . 

If we are interested in just the low-energy, long-wavelength physics, our two-layer system should be describable 
in terms of these operators p"2| . The most general quadratic Hamiltonian we can write down for the system that is 
consistent with the U (1) symmetry in the absence of tunneling takes the form 



H = Jdrjdv' {K {1 (r - r') [S x (r) S x (r') + S y (r) S y (r')] +K±(r- r') S z (r) S z (r')} 



(38) 



The functions Kn and K± are assumed to have a form that will induce a spin density wave, i.e., stripes, in the ground 
state. For example, they could take the (Fourier transformed) form K\\ (q) = p s q 2 , K± (q) = k (— q 2 + <? 4 £ 2 ) where q 
is the wavevector, p s a spin stiffness, and k, S, 2 are positive constants. One can see for n > p s that a uniform spin 
state will be unstable to a state in which the spins tumble spatially; i.e., helimagnetic ordering. The precise form of 
the ground state is difficult to find, even if the spin operators are treated classically; however, qualitatively we know 
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they will have a form similar to the stripe states in our Hartree-Fock analysis. In any case, the results below do not 
depend on any specific choice of K\\ or Kj_, only on the requirement that there is helimagnetic ordering in the ground 
state. 

As is common in a spin- wave analysis [^3, we begin by treating the spins classically. Imposing the constraint 
S x ( r ) 2 + S y (r) 2 + S2 (r) 2 = 1 with a Lagrange multiplier A (r), Eq. (|3§1 ) may be minimized to obtain the three 
equations 

J dx>K\\ (r - r') S x , v (r ; ) - A (r) S x . v (r) 

J dr'K ± (r - r') S z (r') = A (r) S z (r) . (39) 

Eqs. ([39]) together with the constraint equation specify the (classical) ground state. We assume the solutions to these 
equations may be written in the form 

S x (r) = cos 6 (x) = c (x) 
S y (r) = 

S x (r) = sin 0(2;) = s (x) . (40) 
Given the symmetries of H, it is clear that equal energy, inequivalent states can be generated by rotating S in the 



Q 



= 1 



x — y (spin) plane, by translation [0 (x) — ► 8 (x — u) for u a constant], or by rotation [9 (x) — > 6 • rj ,with 

These properties are responsible for the presence of the two gapless modes and their dispersions. 

The spin wave spectrum around this ground state is conveniently found by working in a rotated spin basis, such 
that in the ground state all the spins are aligned along the z axis [E4j. We thus define new spin operators 





s (x) — c(x) 

1 I I S y (r) I . (41) 
c{x) s (x) 



Expanding S' z = J I - S x - S' y 2 w 1 - \ S' x + S' y and making use of Eqs. (gq) , (|g), and fl41f) , to quadratic order 
in S' after some algebra the Hamiltonian may be written as 

H-E = Jdrjdr' [K** (r, r') S' x (r) S' x (r>) + K vv (r - r') S' y (r) S' y (r')] (42) 

where 

K xx ^ r /j = s ^ K \\{r- r') s {x') + c (x) K ± (r - r') c (x') - X (r) 6 (r - r') , (43) 



K yy ( r _ r ') = K{] ( r _ r ') _ A ( r ) § ( r _ r ') _ (44) 

and the ground state energy is 

E = J dv J dr' [c (x) K\\ (r - r') c (a;') 4- s (x) K x (r - r') s (x')] . (45) 

In the classical ground state, S' z = 1. The spin wave approximation amounts to approximating the spin commutation 
relation between S' x (r) and S' (r') by 

[S' x (r) , S'y (r')j =i5(r- r') S' z (r) » iS (r - r') . (46) 

Eq. d42] ) is particularly easy to work with, because these commutation relations allow us to think of S' x as a generalized 
"position" , and S' y as a "momentum" , and in the Hamiltonian they are decoupled. 

An exact computation of the normal modes of Eq. (^) is quite difficult; however, we can understand the basic 
properties of the spectrum through the symmetries of the Hamiltonian and the ground state. Fig. ^ illustrates the 
shapes of s (x) and c (x) in the stripe state for two values oid/i. Taking a to be the distance between stripe centers in 
a single layer (i.e., the width of a full unit cell), it is interesting and important to notice that H is invariant under the 
operation S' (r) — » S' (r + |x); i.e., the primitive unit cell is half the size one expects naively, because s (x) enters the 
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Hamiltonian quadratically, and is invariant under s (x) — ► — s (x + a/2). The discrete translational invariance tells us 
that the normal modes should be expanded in Bloch functions: 

(47) 

(48) 
(49) 

(50) 



a/2 



u n (q) = J tTrh*^ (x) e-^S' x (r) 

Pn (q) = J d 2 rh n , qi (x) e^S' y (r) 
h n ,q x (x + a/2) = h niqa .(x) 

dxh n „ (x) h m _q x (x) — 3mn 



and — 2ix / a < q x < 2n/a defines the effective Brillouin zone. The functions h n>qil . may be chosen so that H — Eq takes 
the form 



\k n (q) u n (-q) u n (q) + - — ^-r^Pn (-q) Pn (q) 
2 2m n (q) 



(51) 



where f2 is the system area, and 

\k n (q) = ^jdrj dr'h* n<qx (x) e"**' if" (r,r') e^'X,,. (x') 

^(qT = hJ dr J dT ' hn ^ {X) e ^ KVV (r " r ° ^"'^.a. 



2m r 



(52) 

From the form of Eq. (|5l|), it is clear that the excitation frequencies of the system are (q) = k n (q) /m n (q) . 
We thus see that there will be gapless (zero) modes whenever k n (q) or l/m„ (q) vanishes. This occurs if there are 
choices for the Bloch functions h n ,q x (x) which satisfy either 



J dr'K xx (r, r') e l 



h n,q m {x') = 



(53) 



J dv'K yv (r 



r^e-^'h* (x') = 



(54) 



Using the symmetries of the ground state, one may find two choices of h, n ^ qm (x) satisfying Eq. 

h* n=Q „-Q (x) =c(x)= cos6» (x) 



or For 



(55) 



it is easily shown that Eq. (|5J) is satisfied. This mode represents a uniform rotation of the ground state spin pattern 
S (r) in the spin x — y plane; i.e., it is associated with the spontaneous phase coherence in the ground state. Making 
use of Eqs. (39), one may also show that Eq. (53) is satisfied for 



--Q/2 (X) 



de 

dx 



-iQ-r/2 



(56) 



for Q = —x. This second zero mode arises due to the translational invariance in H and is a phonon mode. It is 
interesting to note that the two zero modes are found in different parts of the Brillouin zone, the phase mode dispersing 
from the zone center, the phonon from the zone edge of the effective Brillouin zone. In our numerical calculations 
we found both modes dispersing from the zone center. The reason for this is that our numerical technique obliges 
us to work with the naive primitive unit cell, with a resulting Brillouin zone half the size of the one we use in this 
section, so that the phonon mode is folded back to the zone center. It is important to note that since two zero 
modes occur at different wavevectors, they do not mix together and complicate the dispersion uj n (q). This point was 
missed in Ref. p"8| ] , where it was supposed that such mixing would lead to only a single gapless mode with appreciable 
oscillator strength in most response functions. The presence of two gapless modes dispersing from different points in 
the Brillouin zone is precisely what one finds for the spin wave spectrum of an incommensurate helimagnet [ p4| . 
We are left with determining how loq (q) disperses from the two zero modes. In the case of the phase mode, for 



which 1/rao (q) vanishes at q = 0, we may approximate k n —o tqx 
in small powers of q to find 



k 



n=0,g a: =0 



ko 



kg near q = 0, and expand 1/toq (q) 



(57) 
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where 




dr I dv'c (x) K\\ (r - r') (x - x'f c (x') , 
dr I dv'c (x) K\\ (r - r') (y - y'f c (a/) . 



(59) 



(58) 



Thus, the phase mode disperses linearly from q = 0. Notice that if c (x) is only very different from zero in narrow 
regions, as occurs for large layer separations (see Fig. |?]), then c[x) K\\ (r — r')c(x') will be small unless x and x' 

are in the same LCR. Due to the (x — x') 2 factor in pj-, pjr <C pi in this limit; i.e., the dispersion of the phase 
mode perpendicular to the stripes becomes relatively weak. This is precisely the behavior observed in our numerical 
calculations of the collective modes. 

For the phonon mode, it is k n =o,q which vanishes as q — » Q/2. Writing Sq = q — Q/2, it is not difficult to see 
how k n= Q tqx must behave for small \5q\, once one recognizes that the "position" field Mn=o,q=Q/2+<5qi when Fourier 
transformed back to real space represent a displacement perpendicular to the stripes. In this case, the stiffness must 
have the standard smectic form [Q fen=o,q=Q/2+<5q = K±5q^. + Ki>5qy where y is the direction parallel to the stripes, and 
the absence of a Sq y term is a direct result of the rotational symmetry of the Hamiltonian. Kb is a bending modulus 
for the stripes and represents an energy cost for introducing a curvature along them. Writing mo t q=Q/2+Sq ~ m o, we 



Thus, we see the phonon mode disperses linearly with |<5q|, except along the direction parallel to the stripes, for which 
it disperses quadratically. A careful examination of the phonon mode in our numerical results confirms this behavior. 

In closing this section, we note that an observation of the phonon mode dispersion would yield a direct confirmation 
of stripe ordering in this system: the quadratic dispersion along the stripes is indicative of spontaneous smectic 
ordering. 



In a double layer system, the uniform coherent state (UCS) is unstable with respect to the formation of a 
unidirectional and coherent charge density wave state (UCCDW) at a critical value of the interlayer separation 
d — d c (N) .Working in the generalized random-phase approximation (GRPA), we have computed the dispersion rela- 
tions of the low-energy collective modes of the UCCDW in a range of d where this state is expected to be the ground 
state of the 2DEG in the bilayer system. The UCCDW has two Goldstone modes that are respectively related to 
the broken translational symmetry of the stripes and to the broken U(l) symmetry of the coherent state. In the 
long- wavelength limit, the dispersion relations of these modes are consistent with the spin wave dispersion obtained 
in a non-collinear ferromagnet with hclimagnctic ordering. 
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d/l 

FIG. 1. Hartree-Fock energy per electron and coherence as a function of the interlayer separation in the UCS, CWC, UCCDW 
and MUCDW states in Landau level N = 2 and for t = 0. 



0/1=0.7 




FIG. 2. Real space representation of the guiding center and pseudospin densities in the UCCDW for N — 2 and t = 0. 
(a) d/l = 0.7; (b) d/i = 2.0. The width of the linear coherent regions (LCR's) become narrower as the interlayer separation 
increases. 
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FIG. 3. Evolution of the band structure E (X) and of the density of states (DOS) in the UCCDW at N = 2 and t = as a 
function of the interlayer separation. The interstripes separation is approximately a = 6.2£ for these values of d. 
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FIG. 4. Dispersion relations calculated in the GRPA (a) without and (b) with tunneling for the lowest-energy collective 
excitations in the UCCDW in Landau level N = 2 and for d/l = 1.0. These dispersions are for wavevector k in the direction 
of the stripes (k± = 0). The filled diamonds represent the phonon modes (in-phase and out-of-phase) while the lowest two 
branches of the XY mode are represented by the empty squares. 
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FIG. 5. Dispersions of the XY mode of the UCCCW with N = 2 and t = for wavevector (a) parallel (fcj 
perpendicular (fen = 0) to the stripes and for several values of the interlayer separation. 
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FIG. 6. Dispersions of the in-phase phonon mode of the UCCCW with N = 2 and t = for wavevector (a) parallel (fcj 
and (b) perpendicular (fey = 0) to the stripes and for several values of the interlayer separation. 
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FIG. 7. Plots of the functions c(x) and s(x) defined in Eq. ^ for the UCCDW in N = 2 and for t = 0. (a) d/£ = 0.7; (b) 
d/t = 1.5. 
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